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Abstrac t 


Semiempirical  molecular  orbital  calculations  are  reported  in  areas  of 
basic  chemical  Air  Force  research:  explosives,  propellants,  electrochemistry, 
and  the  oxygen-iodine  laser  system.  This  report  is  a  final  report  on  Work 
Unit  2303-F4-03,  "Molecular  Orbital  Calculations  of  Excited  Species".  Work 
will  be  continued  under  a  new  work  unit,  2303-F3-07,  "Theoretical  Chemical 
Calculations  of  Energetic  Materials". 

These  molecular  orbital  calculations  have  given  useful  insights  in  each  of 
these  areas.  In  particular,  we  have  made  substantial  progress  on 
understanding  TNT  thermochemical  decomposition. 
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INTRODUCTION 


Since  its  inception  in  December  1976,  Work  Unit  2303-F4-03,  entitled 
"Molecular  Orbital  Calculations  of  Excited  Species",  has  sought  to  perform 
theoretical  chemical  calculations  on  molecules  of  Air  Force  interest. 
Experimental  chemistry  has  become  quite  expensive,  particularly  in  the  areas 
of  development  of  new  explosives  and  propellants.  In  addition,  some  chemical 
problems,  such  as  characterization  of  transition  states  in  chemical  reactions, 
have  no  experimental  techniques  on  which  they  can  rely  directly  for  answers. 
Our  main  goals  of  these  theoretical  chemical  calculations  are  to  (1)  give 
guidance  to  experiments  designed  to  develop  new  energetic  molecules,  and  (2) 
explain  experimental  evidence  available  in  terms  of  possible  chemical  reaction 
mechanisms. 

Work  has  been  done  in  the  seven  years  of  this  work  unit  in  a  variety  of 
areas.  This  final  report  is  meant  to  be  a  summation  of  all  previously 
published  work  from  the  work  unit  as  well  as  a  report  on  some  as  yet 
unreported  calculations,  those  being  in  the  explosives  decomposition  area. 

The  report  is  broken  into  sections  according  to  the  chemical  system  being 
studied. 

In  addition  to  the  work  reported  here,  there  has  been  other  work  in  the 
theoretical  chemistry  area  supporting  other  tasks.  This  other  work  includes  a 

study  of  the 'potential  of  certain  energetic  plasticizer  molecules  in 

1  .  .  .2 

propellants  and  the  decomposition  modes  of  tetrazines. 

There  is  much  left  to  be  done,  again',  particularly  in  the  explosives 
decomposition  area.  This  work  unit  (2303-F4-03)  is  being  closed  out  because 
of  the  termination  of  task  2303-F4  (Molecular  Dynamics).  Because  most  ot  the 
calculations  have  been  done  on  energetic  :.ater:  1 1  s  r  new  wo  unit 
(2303-F3-07) ,  "Theoretical  Calculations  of  Energetic  Materials",  is  being 
initiated  in  the  Energetic  Materials  Task  (2303-F3)  in  FY  84.  The  same  type 
of  theoretical  chemical  calculations  will  be  continued  in  that  work  unit. 
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Chapter  2 

THEORETICAL  METHODS 

All  of  the  work  performed  in  this  work  unit  has  been  with  the  semiempirica  1 
methods  developed  by  Professor  M.  J.  S.  Dewar  and  coworkers  at  the  University 
of  Texas  in  Austin,  TX.  Just  recently  we  have  begun  some  preliminary  work 
with  the  jib  initio  program  GAUSSIAN  80,  but  that  work  will  be  reported  in  the 
final  report  of  University  Resident  Research  Program  Fellow  Dr.  Almon  G. 
Turner. 

3 

The  first  method  utilized  was  the  MINDO/3  scheme  ,  It  is  a 
semiempirical  molecular  orbital  program  based  on  an  Intermediate  Neglect  of 
Differential  Overlap  (INDO),  approximation  and  has  the  capability  of 

optimizing  geometries  upon  request.  Dewar's  more  recent  program,  called 

4  .  .  , 

MNDO  ,  is  a  semiempirical  method  based  on  a  Neglect  of  Diatomic  Differential 
Overlap  (NDDO)  approximation ,  and  it,  like  MINDO/3  can  optimize  molecule 
geometries  upon  request. 

Each  method  calculates  a  variety  of  molecular  properties  of  interest, 
including  equilibrium  molecular  geometries,  heats  of  formation,  dipole 
moments,  atomic  charges,  and  gradients  along  the  potential  surface.  In 
general,  MNDO  results  have  proved  superior  to  MINDO/3,  except  for  certain 
classes  of  molecules,  the  most  notable  of  which  is  nitro-containing 
compounds.  Since  this  class  of  compounds  is  very  important  in  our  explosives 
work,  we  thought  it  wise  to  do  calculations  with  both  methods  and  compare 
results. 

Related  programs  also  allow  searches  for  transition  states  along  a 
reaction  pathway,  r^rce  constant  calculations,  and  thermodynamic 
calculat  ions.  The  force  constant  calculations  are  important  in  showing  where 
one  is  on  a  molecular  potential  surface.  At  a  stationary  point,  all  gradients 
are  very  close  to  zero.  If  the  stationary  point  is  a  stable  molecule,  all  of 
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the  3N  -  6  force  constants  corresponding  to  vibrational  motion  will  be 
positive.  If,  however,  the  stationary  point  has  one,  and  only  one,  negative 
force  constant,  it  must  be  a  transition  state,  with  the  negative  force 
constant  indicating  that  that  mode  of  motion  is  the  reaction  coordinate.  Two 
or  more  negative  force  constants  indicate  a  "hill"  on  the  surface,  which 
cannot  be  a  transition  state.  The  results  of  the  force  constant  calculations 
can  be  input  into  a  thermodynamic  program  in  order  to  get  partition  functions, 
entropies,  and  heat  capacities. 

Within  each  method  (MNDO  and  MINDO/3),  it  is  sometimes  advisable  to  use  a 
Unrestricted  Hartree-Fock  (UHF)  technique  if  the  system  of  interest  is  a 
reaction  path  or  transition  state  in  which  radicals  are  produced  from  closed 
shell  molecules.  Also,  UHF  optimizations  of  transition  states  are  often 
faster,  and  the  optimized  geometry  of  the  UHF  calculation  is  often  used  to 
begin  the  RHF  or  half-electron  calculation. 

Recently,  Stewart  and  Dewar  have  combined  all  of  these  options  into  one 
Molecular  Oirbital  Package  (MOPAC)^.  The  combined  package  has  greatly 
improved  ease  of  use  and  is  quite  powerful  in  its  capabilities. 
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ELECTROCHEMICAL  CALCULATIONS 


One  important  measurable  quantity  in  electrochemical  behavior  of  solutes 
is  the  half-wave  potential.  The  quantity  is  a  measure  of  how  difficult  it  is 
to  remove  an  electron  (or  oxidize)  the  species  in  solution.  Since  the 
calculated  molecular  orbital  energies  are  directly  related  to  the  molecular 
ionization  energy  via  Koopmans'  theorem^,  we  sought  to  use  the  calculated 
molecular  orbital  energies  from  MNDO  to  predict  the  oxidation  half-wave 
potentials  in  a  room-temperature  organic  molten  salt  melt.^ 

We  found  good  correlation  in  the  study  between  the  measured  oxidation 
half-wave  potentials  and  the  calculated  conjation  energies  for  a  series  of 
aromatic  compounds,  including  some  nitroaromatics  in  acetonitrile.  The 
half-wave  potential  could  be  calculated  from  the  predicted  ionization  energy 
according  to: 

E1/2  =  0.725  Ej  -  4.95 

Ej^2  is  the  measured  half-wave  potential  in  volts  and  E^  is  the  calculated 
ionization  potential  in  ev. 

In  the  room  temperature  molten  salt,  we  found  a  similar  correlation  for 
the  same  set  of  compounds: 

E  =  0.725  ET  -  4.72 
m  I 

E  is  the  half-wave  potential  in  the  melt  in  volts, 
m 

The  prediction  for  thianthrene  was  not  as  good  as  the  other  species,  being 
off  from  the  experimented  value  by  about  1  volt.  This  discrepenoy  may  be  due 
to  the  fact  that  the  highest  occupied  molecular  orbital  (HOMO)  for  thianthrene 
is  a  lone-pair  orbital  on  the  sulfur,  while  for  all  of  the  other  molecules, 
the  HOMO  was  an  aromatic  Tt-type  orbital. 

The  results  have  been  used  to  predict  half-wave  potentials  for  other 
compounds  to  determine  their  suitability  as  electrochemical  species  in  the 


room  temperature  molten  salt.  In  some  cases,  the  molecule  is  spontaneously 
oxidized  by  the  melt.  According  to  the  theoretical  correlation,  this  would 
happen  whenever  the  predicted  ionization  potential  is  less  than  6.51  ev , 
resulting  in  a  negative  half-wave  potential. 

The  other  end  of  the  electrochemical  window  is  al6o  important.  If  the 
species  is  spontaneously  reduced  by  the  melt,  it  will  also  be  unstable.  The 
reduction  potential  of  the  molecule  is  related  to  its  electron  affinity,  its 
desire  to  acquire  an  electron.  In  theoretical  calculations,  this  electron 
affinity  shows  up  as  the  energy  of  the  lowest  unoccupied  molecular  orbital 
(LUMO).  We  did  some  calculations  to  confirm  the  usefulness  of  MNDO  to 
accurately  predict  the  electron  affinities  of  a  group  of  molecules,  and  found 

g 

good  agreement.  These  correlations  have  proven  quite  useful  in 
establishing  the  electrochemical  window  for  a  variety  of  melt  components  and 
species  on  solution. 
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Chapter  4 

Aluminum  Calculations 

The  Air  Force  is  interested  in  chemical  properties  of  rocket  exhaust  plume 
species  in  order  to  determine  the  plume's  radiation,  contamination,  and  other 
effects.  Many  current  rocket  engines  use  fuels  co>  ....ining  aluminum. 

Therefore,  in  order  to  calculate  properties  of  various  aluminum-containing 
molecules,  it  was  necessary  to  determine  MNDO  parameters  r  r  aluminum  and 
assess  the  accuracy  of  MNDO  calculations  for  these  compounds. 

The  parameterization  work  was  done  by  Henry  Rzepa  at  the  University  of 

Texas  (Dewar's  Croup)  and  the  assessment  by  Davis,  Guidry,  and  Williams  here 

9 

at  FJSRL.  Both  items  were  combined  into  a  single  journal  article.  The 
article  also  included  some  work  on  boron-containing  compounds. 

The  results  indicated  that  the  results  for  the  aluminum-containing 
molecules  were  clearly  in  excess  of  those  for  compounds  containing  only  C,  H, 

0,  and  N.  The  average  error  in  the  heats  of  formation  for  the  aluminum- 
containing  molecules  was  20.3  kcal  mol  Several  trends  were  noted, 

however.  The  largest  errors  resulted  from  molecules  which  contained  oxygen  as 
well  as  the  aluminum  (average  error  of  32.6  kcal  mol  In  addition, 

radicals  were  predicted  to  he  particularly  poor  (average  error  of  31,2  kcal 
m°l  )•  As  the  molecule  became  larger,  the  predictions  became  progressively 
too  endothermic. 

Geometry  predictions  fared  reasonably  well.  The  only  obvious  trend  was 
the  under pred ic t i on  of  the  aluminum  bond  lengths.  The  general  shapes  of 
molecules  were  well  predicted.  The  errors  in  geometry  indicate  that  d 
orbitals  may  be  required  to  accurately  predict  the  geometry  ot 
aluminum-containing  molecules. 

Thus  the  results  for  aluminum,  while  they  do  need  some  improvement,  are 
still  accurate  enough  to  provide  guidance  in  explaining  chemical  reaction 
mechanisms  and  calculating  properties  of  aluminum-containing  molecules. 
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Oxygen  Singlet-Delta  Calculat ions 


One  potential  new  chemical  laser  system  in  the  near-IR  involves  reaction 

of  basic  hydrogen  peroxide  with  gaseous  chlorine  to  produce  the 

electronically-excited  02(^6).  The  02^4)  can  pump  an  atomic 

iodine  laser^'  One  of  the  major  unknowns  in  developing  these  reactions 

into  an  operational  chemical  laser  is  the  mechanism  of  the  reaction  in  the 

chemical  generator  to  produce  the  OjC^A).  We  have  used  MNDO 

12 

calculations  to  suggest  plausible  mechanisms. 

The  overall  reaction  is 

Cl2(g)  +  H202(aq)  ♦  20H~(aq)  *  02<1A)(g)  +  21^0  2Cl~(aq) 

The  efficiency  of  the  O2  generator  is  very  pH-dependent ,  decreasing  as  the 
pH  decreases.  In  addition,  the  system  is  complicated  considerably  by  the  two- 
phase  nature  of  the  process. 

Specifically,  our  goals  with  MNDO  were  (1)  to  test  the  feasibility  of 

proposed  reaction  intermediates  (2)  to  determine  whether  energet ical ly- 

feasible  intermediates  will  decompose  into  02(^A)  and  other  products 

and,  (3)  to  assess  the  likelihood  of  intersystem  crossing  to  the  triplet 

3 

surface  to  produce  02(  E) . 

First  we  considered  possible  species  present  in  the  solution  before  the 
Cljig)  was  added.  Several  acid-base  equilibria  are  present,  but  both  the 
calculations  and  available  experimental  data  indicate  that  the  predominant 
species  in  the  basic  hydrogen  peroxide  solution  is  the  hydroperoxide  ion, 

00H  .  Of  lesser  importance  is  the  hydroxide  ion,  OH  ,  and  some  0C1  may 
be  present  after  the  Clj  is  introduced  to  the  solution. 

Now  we  considered  reactions  of  the  C^Cg)  with  these  solution  species. 

Clj  reacts  with  00H  without  activation  (according  to  the  calculation) 
either  to  directly  produce  Oj^A),  HC1,  and  Cl  ,  or  to  form  a  complex 


HOOClj  •  The  direct  reaction  occurs  only  if  the  attack  of  the  Cl2  is 
constrained  to  occur  only  on  the  H  of  the  00H  ;  otherwise  the  complex 
forms.  The  complex  is  predicted  to  be  a  stable  intermediate,  but  there  is 
only  a  fairly  small  barrier  to  dissociation  of  a  chloride  ion  to  form  H00C1  + 
Cl". 

In  either  case  (H00C1  or  H00C12  ),  attack  by  any  basic  species  in 
solution  is  predicted  to  break  up  the  complex  and  form  OjC^A)  without 
activation.  Thus  (^(^A)  can  be  formed  without  activation  either  by 
direct  reaction  of  Cl2  with  00H  or  through  the  formation  of  an 
intermediate  complex. 

Finally,  we  considered  the  possibility  of  intersystem  crossing  to  the 
triplet  state,  since  this  would  decrease  the  yield  of  (^(^A).  The 
calculations  producted  no  feasible  path  for  intersystem  crossing  in  basic 
solution.  In  acidic  solution,  some  H0C1  may  form,  and  there  is  a  possibility 
of  intersystem  crossing  in  the  reaction  of  OOH  and  H0C1 . 

The  results  showed  the  feasibility  of  using  MNDO  to  investigate  chemical 
reaction  mechanisms  in  solution. 


8 


r 


Chapter  6 

Decomposition  of  Explosives  and  Propellants 


Introduct ion 

Most  of  the  work  in  this  area  has  focussed  on  the  thermochemical 
decomposition  of  TNT,  TNT-model  compounds,  RDX,  and  RDX-model  compounds. 
However,  we  also  studied  the  decomposition  modes  of  hydrazine  (^H^)  in 
spacecraft  fuel  tanks  at  the  request  of  AFWL.  We  will  first  report  on  these 
calculations. 

Decomposition  of  N- 

AFWL  needed  activation  energies  for  all  possible  decomposition  modes  of 
We  studied  breaking  either  the  N-N  bond  or  the  N-H  bond  for 
hydrazine,  the  hydrazine  cation,  and  the  hydrazine  dication.  Results 
indicated  that  removal  of  a  hydrogen  atom  had  an  activation  energy  of  77.6 
kcal  mol  *  for  hydrazine,  75.3  kcal  mol  1  for  the  cation,  and  48.1  kcal 
mol  *  for  the  dication.  Breaking  the  N-N  bond  had  an  activation  energy  of 
68.0  kcal  mol  *  for  hydrazine,  127.2  kcal  mol  ^  for  the  cation,  and  108.3 
kcal  mol  ^  for  the  dication.  These  results  were  transmitted  to  AFWL  and 
used  in  their  modelling  studies. 
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TNT  Decomposition 

For  several  years  the  kinetics  and  mechanism  of  the  thermochemical 

13 

decomposition  of  TNT  has  been  studied  under  Work  Unit  2303-F3-05.  We  have 

14 

sought  to  theoretically  .ivestigate  this  topic  as  well.  Because  of  the 
size  of  the  TNT  molecule,  and  hence  its  long  computation  times,  we  have  used 
1-nitropropene  as  a  model  compound  for  TNT: 


TNT 


1-NITROPROPENE 


Note  that  the  cis-l-nitropropene  has  the  same  structural  features,  if  one  is 
investigating  reactions  which  involve  the  methyl  and  nitro  groups,  as  TNT. 

The  TNT  decomposition  mechanism  has  been  discussed  for  years  in  the 
literature.  Based  on  these  discussions,  we  chose  the  following  six  reactions 
as  possible  first  steps  for  the  TNT  decomposition: 
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(1)  [r* ramolecular  Hydrogen  Transfer 


CH, 


l-nitropropene  ♦ 


2\  /  N°zH 

/•c. 

K  H 


(2)  Intermolecular  Hydrogen  Transfer 


2  (l-nitropropene) 


NO. 


CH2v  /,W2 
C  -  c 
/  \ 

K 


H 


(3)  Intramolecular  Oxygen  Insertion 


l-nitropropene 


hoch  n-o 

/ 


c  -  c 
✓  \ 

H  H 


(4)  Intermolecular  Oxygen  Insertion 


2  (l-nitropropene) 


CH 


NO 

3\  ' 

/' c- 

H  » 


(5)  Nitro-Nitrite  Rearrangement 


l-nitropropene 


CH 


3v  / 
c  -  c 


O-N-O 


(6)  Nitro  Group  Removal 


l-nitropropene  * 


Clt 


3s 

/C-C; 

it  >i 


CH2  /°2U 

m  C 

/ 


H 


CH.OH  NO 

‘S  *  1 


c  -  c 

/  V 
H  H 


♦  NO, 
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If  I-nitropropene  is  to  be  a  good  model  for  TNT,  we  expect  the  AH0,s  for 
these  reactions  to  be  similar  for  the  two  systems.  Since  calculating  heats  of 
formation  for  known  stationary  points  (reactants  and  products)  is  not  overly 
time-consuming,  we  calculated  AH°'s  for  all  of  these  six  reactions  both  with 
TNT  and  1-nitropropene.  The  results  are  given  in  Table  I. 
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TABLE  X.  STANDARD  ENTHALPIES  OF  REACTION  (kcal  mol'1) 


Reaction 


MNDO 


MINDO/3 


1-Nitropropene 

TNT 

A. 

1-Nitropropene 

TNT 

Intra  H 
Transfer 

8.6 

16.6 

8.0 

7.6 

19.8 

12.2 

Inter  H 
Transfer 

16.2 

16.8 

0.6 

.38.2 

38.7 

0.5 

Intra  0 
Insertion 

-49.5 

-49.9 

-0.4 

-18.7 

1 

H- * 

00 

0.3 

Inter  0 
Insertion 

-51.3 

-52.4 

-1.1 

-16.1 

-21.9 

-5.8 

Nitro-Nitrite 

Rearrange¬ 

ment 

-38.9 

-41.1 

-2.2 

-11.8 

-7.5 

4.3 

Nitro  Group 
Removal 

29.8 

26.9 

-2.9 

41.4 

62.1 

20.  7 

13 


We  note  that  in  general  there  is  good  agreement  between  the 
1-nitropropene  result  and  the  TNT  result  for  either  method.  There  is 
a  large  discrepancy,  however,  between  the  MNDO  and  MINDO/3  results  for 
most  reactions.  When  they  differ,  we  would  expect  the  MINDO/3  results 
to  be  better,  since  MNDO  gives  poor  results  for  heats  of  formation  for 
nitro-containing  compounds. ^  The  only  really  bad  result  that 
MINDO/3  gives  is  for  NC^,  making  it  about  25  kcal  mol  ^  too 
Stable.  Thus  the  AH°  for  the  nitro  group  removal  will  be 
underestimated.  Since  MNDO  generally  predicts  nitro  compounds  too 
unstable,  it  will  in  general  make  the  predicted  AH0  too  exothermic 
for  any  reaction  which  involves  changing  the  nitro  group  in  some  way. 
This  is  borne  out  by  comparing  the  MNDO  and  MINDO/3  results. 

Since  1-nitropropene  does  indeed  correlate  well  with  TNT  in  terms 
of  AH®'s,  it  should  also  serve  as  a  good  model  for  the  activation 
enthalpies  of  these  processes.  Transition  states  were  located  by 
defining  a  reaction  coordinate  for  each  reaction,  finding  the  point  of 
maximum  enthalpy,  and  then  using  a  transition  state  location  program 
which  minimizes  the  sum  of  squares  of  the  gradients.  Each  transition 
state  was  identified  as  a  transition  state  by  running  a  force  constant 
calculation  which  produced  one,  and  only  one,  negative  force 
constant.  Table  II  gives  the  calculated  activation  enthalpies  for 
each  reaction. 
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TABLE  II.  ACTIVATION  ENTHALPIES  (kcal  mol"1)  FOR  I-NITROPROPENE  SYSTEM 


Reaction 


MNDO 


MINDO/3 


Intra  H 

Transfer  63.7  (59.3)  43.4  (43. 4^ 


Inter  H 

Transfer  —  (50.4)  —  (39.6) 


Intra  0 

Insertion  127.7  49.9 


Inter  0 

Insertion  108.7  44.4 


Nitro-Nitrite 
Rearrange¬ 
ment  51.6  35.0 


Nitro  Group* 

Removal  29.8  41.4 


(  )  -  UHF  Calculation 

*  No  true  transition  state;  6H°  assumed  to  be  the  reaction  AH°. 
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The  transition  state  obtained  for  the  intramolecular  hydrogen  atom 


transfer  of  a  hydrogen  atom  from  the  CH^  end  of  the  three  carbon  chain  to 
the  nitro  group  is  a  cyclic  structure  wherein  the  hydrogen  atom-carbon  atom 
bond  has  been  elongated  from  1.10A  to  1.34A  (bond  order  of  0.42).  The 
initial  stages  of  O-H  bond  formation  can  be  identified  with  the  hydrogen 
oxygen  distance  of  1.17A  (bond  order  of  0.44).  The  transition  state  is 
non-planar  with  the  hydrogen  being  transferred  about  0.5A  out  of  the  plane 
formed  by  the  three  carbon  atoms.  The  force  constant  calculation  gave 
excellent  results:  one  large  negative  eigenvalue  for  the  reaction  coordinate 
and  six  close-to-zero  results  for  the  translations  and  rotations  (all  less 
than  6  cm  ^).  The  reaction  coordinate  corresponded  primarily  to  motion  of 
the  hydrogen  being  transferred. 

The  calculated  activation  enthalpy  was  somewhat  higher  for  MNDO  than 
MINDO/3  (63.7  kcal  mol  ^  compared  to  43.4  kcal  mol  ^).  This  is  consistent 
with  MNDO  normally  giving  slightly  high  activation  enthalpies  for  hydrogen 
transfer  reactions.  The  MINDO/3  result  is  consistent  with  the  experimental 
activation  energy  and  thus  this  process  cannot  be  ruled  out  as  a  first  step  in 
the  thermochemical  decomposition  on  the  basis  of  this  calculated  activation 
energy.  We  would  expect  the  TNT  activation  enthalpy  to  be  several  kcal  mol 
higher;  however,  because  of  the  loss  of  the  aromatic  character  of  the  ring  as 
TNT  converts  to  the  product  state.  This  effect,  of  course,  will  not  show  up 
in  our  1-nitropropene  model  compound.  Because  of  the  possible  diradical 
character  of  the  product,  we  also  optimized  the  transition  state  using  an  UHF 
type  calculation.  For  MINDO/3,  no  spin  imbalance  was  calculated  and  the 
results  were  the  same  as  the  RHF  results.  MNDO  gave  a  slightly  lower  value 
for  the  activation  enthalpy  with  UHF,  but  the  difference  was  not  significant. 

The  transition  state  for  the  intermolecular  hydrogen  transfer  was 
calculated  using  UHF.  The  product  of  this  reaction  is  two  radicals  and 
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attempts  to  optimize  an  RHF  transition  state  resulted,  as  could  be  expected, 
in  the  separated  reactants.  In  contrast  to  the  intramolecular  hydrogen 
transfer  transition  state,  which  was  almost  symmetrical  in  terms  of  the  C-H 
and  0-H  bond  orders,  this  intermolecular  transition  state  shows  that  the  C-H 
bond  is  more  than  half-broken  (C-H  distance  is  1.36A,  with  a  bond  order  of 
0.30)  while  the  O-H  bond  is  more  appreciably  formed  (O-H  distance  is  1 . 0 7 A , 
with  a  bond  order  of  0.53).  The  C-H-0  bond  angle  is  166°,  approaching 
linearity.  The  force  constant  calculation  gave  one  large  negative  eigenvalue 
for  the  reaction  coordinate,  but  the  translations  and  vibrations  did  mix 
appreciably  with  the  low-lying  vibrations  (frequencies  for  the  translations 
and  rotations  exceeded  100  cm  ^).  This  mixing  is  not  surprising  considering 
the  loosely  bound  nature  of  this  bimolecular  transition  state,  and  it  should 
not  qualitatively  affect  the  conclusions.  Once  again,  the  reaction  coordinate 
indicated  the  motion  to  consist  primarily  of  hydrogen  atom  translations. 

Again,  the  calculated  activation  enthalpy  for  MND0  is  somewhat  higher  than 
for  MINDO/3  (50.4  kcal  mol  vs.  39.6  kcal  mol  ^),  following  the  trend 
that  MNDO  overpredicts  activation  enthalpies  lor  hydrogen  transfer  reactions. 
The  MINDO/3  result  for  this  reaction  is  also  consistent  with  the  experimental 
activation  energy,  although  MINDO/3  UHF  results  typically  may  underestimate 
the  activation  enthalpy  by  several  kcal  mol  Even  so,  the  agreement  with 

the  experimental  number  is  quite  good. 

In  the  MINDO/3  transition  state  for  the  intramolecular  oxygen  insertion 
from  the  nitro  group  into  a  C-H  bond  of  the  methyl  group,  the  N-0  bond  is 
stretched  (distance  of  1.35A,  bond  order  of  0.66),  tiie  C-H  bond  is  stretched 
(distance  of  1.25A,  bond  order  of  0.59),  thus  0-0  bond  is  paiiiallv  formed 
(distance  of  l . 6  3 A ,  bond  order  of  0.29),  and  the  O-H  bond  is  partially 
formed  (distance  of  L.22A,,bond  order  of  0.29).  The  geometry  of  the 
transition  state  is  slightly  non-planai ,  with  the  inserted  oxygen  atom  and  the 
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hydrogen  atom  both  about  O.lA  out  of  the  plane  formed  by  the  three  carbon 
atoms.  Indeed,  the  first  attempt  to  optimize  this  transition  state  with 
MINDO/3  produced  a  planar  geometry,  but  the  presence  of  more  than  one  negative 
force  constant  indicated  that  the  structure  was  not  a  true  transition  state. 
The  force  constant  calculation  on  the  true  transition  state  had  only  one 
negative  force  constant,  with  each  of  the  translational  and  rotational  modes 
having  frequencies  less  than  10  cm  The  reaction  coordinate  corresponded 
to  motion  of  all  of  the  atoms  involved  in  either  bond-making  or  bond-breaking, 
with  the  oxygen  and  the  hydrogen  showing  the  largest  effects.  Interestingly, 
there  was  virtually  no  difference  in  the  calculated  heats  of  formation  between 
the  planar  and  non-planar  structure.  The  MNDO  result  was  not  quite  as  good. 

In  spite  of  repeated  attempts  to  displace  the  stationary  point  calculated  by 
MNDO  from  planarity,  we  were  not  successful.  The  result  is  a  planar 
stationary  point  which  we  believe  to  be  very  near  the  true  transition  state, 
but  which  nevertheless  has  two  negative  eigenvalues  in  the  force  constant 
matrix:  one  (corresponding  to  the  true  reaction  coordinate)  of  1733  cm  ^  and 
the  other  (corresponding  primarily  to  torsional  motion  of  the  hydrogen  in  and 
out  of  the  plane  of  the  molecule)  of  303  cm 

There  is  a  large  difference  between  the  calculated  enthalpies  of 
activation  for  MNDO  and  MINDO/3  (127.7  kcal  mol  ^  vs.  49.9  kcal  mol  ^), 
even  though  the  overall  geometries  are  similar.  Based  on  previously-discussed 
rationale,  the  MINDO/3  results  are  more  believable.  The  MINDO/3  enthalpy  of 
activation  is  slightly  larger  than  the  experimented  activation  energy,  but  the 
10  kcal  mol  ^  difference  does  not  completely  rule  this  step  out  as  a 
possible  first  step.  We  think,  however,  that  it  is  less  likely  than  the 
hydrogen  transfer  steps  previously  discussed. 

Compared  with  the  intramolecular  insertion  transition  state,  in  the  case 
of  the  intermolecular  oxygen  insertion,  the  reaction  has  proceeded  more  toward 
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product;  the  N-0  bond  is  stretched  (distance  of  1.38A,  bond  order  of  0.75)  , 
the  C-H  bond  is  stretched  (distance  of  1.29A,  bond  order  of  0.98),  the  C-0 
bond  is  partially  formed  (distance  of  1.66A,  bond  order  of  0.35),  and  the 
0-H  bond  is  partially  formed  (distance  of  1 . 1  lA  ,  bond  order  of  0.39).  The 
N-O-C  bond  angle  is  158°,  approaching  linearity.  The  force  constant 
calculation  gives  a  single  negative  force  constant,  with  rather  poor 
separation  erf  the  translational  and  rotational  modes  (maximum  frequency  of  69 
cm  ^).  Again,  this  mixing  with  low-lying  vibrations  is  not  uncommon  for 
bimolecular  transition  states.  The  reaction  coordinate  consists  of  motion 
primarily  by  the  atoms  involved  in  bond-breaking  and  bond-making,  but  there 
are  minor  contributions  from  many  of  the  other  modes. 

Like  the  intramolecular  oxygen  insertion,  there  is  a  large  discrepancy 
between  MNDO  and  MINDO/3  in  the  calculated  enthalpies  of  activation  for  this 
intermolecular  oxygen  insertion  (108.7  kcal  mol  ^  vs.  44.4  kcal  mol  *). 

Again,  assuming  the  MINDO/3  value  to  be  more  accurate,  this  calculated 
activation  enthalpy  is  consistent  with  the  experimental  activation  energy. 

Thi^  process  is  quite  likely  to  be  taking  place  in  the  initial  stages  of  TNT 
decomposition. 

In  the  nitro-nitrite  rearrangement,  one  of  the  oxygens  inserts  into  the 
C-N  bond,  resulting  in  the  structure  C-0 -N-O.  The  transition  state  shows 
stretching  of  the  C-N  bond  (distance  of  1.45A,  bond  order  of  0.71), 
stretching  of  the  N-0  partial  double  bond  in  the  nitro  group  toward  the  single 
bond  in  the  product  (distance  of  1.30A,  bond  order  of  1.19),  and  partial 
formation  of  the  C-0  bond  (distance  of  1.53A,  bond  order  of  0.45).  The 
force  constant  calculation  gives  a  single  negative  eigenvalue  ot  frequency  '05 
cm  with  good  separation  of  the  translational  and  rotational  modes 
(maximum  frequency  of  9  cm  ^).  The  reaction  coordinate  consists  of  movement 
primarily  by  the  three  atoms  involved,  C,  N,  and  0. 
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The  results  for  the  activation  enthalpies  are  35.0  kcal  mol  *  and  51.6 
kcal  mol  for  MINDO/3  and  MNDO,  respectively.  Both  values  are  not  far 
removed  from  the  experimental  activation  energy,  making  this  process  a 
feasible  first  step  on  an  activation  enthalpy  basis. 

The  last  process  considered  was  a  simple  C-N  bond  fission,  resulting  in 
the  C  H,.  radical  and  N0^ .  A  reaction  path  calculation  failed  to  produce 
a  transition  state,  meaning  that  the  reverse  reaction  should  process  without 
activation,  a  common  occurrance  for  a  radical  recombination  reaction.  Thus, 
the  activation  enthalpy  and  the  reaction  enthalpy  are  one  and  the  same.  The 
results  are  29.8  kcal  mol  1  for  MNDO  and  41.4  kcal  mol  1  for  MINDO/3. 

Both  values  make  the  C-N  bond  scission  look  feasible  when  compared  with  the 
experimental  activation  energy,  but  we  feel  that  both  numbers  are  in  error. 
MNDO  results  for  heats  of  formation  of  nitroaromatics  generally  are  far  too 
endothermic,1^  destabilizing  the  reactant  and  making  the  activation  enthalpy 
too  small.  On  the  other  hand,  MINDO/3  gravely  overstabilizes  NOj  by  about 
25  kcal  mol  1 ,  again  making  the  activation  enthalpy  too  small.  If  one 
corrects  the  MINDO/3  result  for  this  known  error  for  NO^,  an  activation 
enthalpy  of  about  65  kcal  mol  1  is  obtained.  This  number  agrees  with  the  70 
kcal  mol  value  one  obtains  for  the  C-N  bond  strength  in  nitrobenzene  if 
the  experimental  heats  of  formation  for  nitrobenzene  and  N0„,  and  the 
MNDO-calculated  value  for  phenyl  radical  are  used.  This  70  kcal  mol  1  value 
is  also  in  the  area  of  a  C-N  bond  strength  according  to  conventional  organic 
chemistry  wisdom.  Thus,  we  feel  that  a  simple  C-N  bond  fission  is  not 
possible  as  a  first  step  in  the  thermochemical  decomposition  of  TNT. 

Cis-l-nitropropene  is  a  good  theoretical  model  for  studying  TNT 
thermochemical  decomposition.  It  reproduces  all  of  the  major  structural 
features  of  a  nitro  group  being  ortho  to  the  methyl  group  in  TNT.  Of  all  of 
the  possible  reaction  steps  we  considered,  only  the  intramolecular  hydrogen 


transfer  changes  the  aromaticity  of  the  TNT  system,  making  our  1-ni tropropene 
model  not  quite  as  good  for  that  case.  Even  so,  we  are  able  to  deduce 
qualitative  trends  of  this  effect. 

MINDO/3  results  do  appear  to  be  more  reliable  than  the  MNDO  results.  We 
cannot  explain  the  abnormally  high  values  of  the  activation  enthalpies  that 
MNDO  gives  for  the  oxygen  transfer  steps.  It  is  quite  feasible  that,  in  spite 
of  repeated  attempts,  we  have  still  not  obtained  the  transition  states  for 
these  processes  with  MNDO.  In  the  other  cases,  the  MINDO/3  and  MNDO  results 
are  similar. 

Our  original  goal  was  to  draw  conclusions  about  possible  first  steps  of 
TNT  thermochemical  decomposition  by  comparing  our  theoretical  results  with  the 
EPR  and  IDSC  results.  On  an  activation  energy  basis,  several  processes  appear 
to  be  likely:  both  hydrogen  transfers,  both  oxygen  transfers,  and  the 
nitro-nitrite  rearrangement.  The  fairly  large  isotope  effect,  however,  is 
difficult  to  explain  with  anything  but  the  two  hydrogen  transfer  steps, 
although  we  are  currently  calculating  deuterium  isotope  effects  for  all  of 
these  steps.  The  EPR  signals  indicate  that  radicals  must  be  involved,  and  the 
spectrum  shows  two  nitrogens  and  five  hydrogens.  The  intermolecular  hydrogen 
transfer  has  the  attraction  that  it  produces  radicals  directly,  but  not  with 
the  required  two  nitrogens.  MONO  loss  from  the  hydrogen  transfer  product  does 
result  in  a  radical  with  the  required  two  nitrogens  and  five  hydrogens,  and 
further  work  is  in  progress  to  determine  the  feasibility  of  these  follow-on 
steps.  On  the  other  hand,  intramolecular  hydrogen  transfer  followed  by  HONO 
loss  would  result  in  a  radical,  but  with  two  nitrogens  and  only  four  hydrogens. 

The  intermolecular  hydrogen  transfer  process  is  consistent  with  all  of  the 
known  experimental  facts  for  TNT  thermochemical  decomposition.  Thus,  we  feel 
that  it  is  a  major  process  early  in  the  decomposition.  The  low  activation 
enthalpies  calculated  for  several  other  steps,  however,  make  it  probable  that 
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these  other  processes  go  on  just  as  well.  Indeed,  all  of  previous  studies  on 
TNT  thermochemical  decomposition  indicate  that  the  system  is  quite 
complicated,  particularly  in  light  of  the  multitude  of  intermediate  products. 
Our  calculations  certainly  indicate  the  complexity  of  the  system. 
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Decomposition  of  RDX 

The  chemistry  of  the  heterocyclic  nitramine  1 , 3 , 5-t r ini t ro-1 , 3 , 5- 
tr iazacyc lohexane  (RDX)  is  of  great  interest  because  of  its  use  as  an 
explosive. 


NO-, 


N 

✓  \ 


CHi 


NO-, 


CHo 

I 

CH2  N02 


RDX 


no2 

I 

N 

✓  \ 

ch3  ch3 


DMNA 


The  purpose  of  this  work  was  to  perform  computer  modeling  of  possible 
thermal  decomposition  paths  of  dimethy lnitramine  (DMNA).  The  reason  that 
DMNA  was  chosen  was  because  of  its  similarity  to  RDX  while  being  much 
simpler  to  model  with  the  computer.  The  following  list  explains  all  the 
decomposition  paths  which  have  been  worked  on  and  which  have  been 
considered  of  interest  to  this  date. 

1.  Breaking  the  N-N  bond  and  dissociating  nitrogen  dioxide. 

2.  Breaking  the  N-N  bond  of  the  anion  and  dissociating  NO,,  . 

3.  Breaking  a  C-N  bond  and  dissociating  a  CH^  group. 

A.  Breaking  a  N-0  and  dissociating  an  oxygen  atom. 

5.  Breaking  a  C-H  bond  and  dissociating  a  hydrogen  atom. 

6.  The  addition  of  a  hydrogen  atom  to  the  NO2  group. 

7.  The  intermolecular  transfer  of  a  hydrogen  atom  between  2  DMNA 

molecules  creating  an  0-H  bond  and  breaking  a  C-H  bond. 

8.  Breaking  the  N-N  bond  and  releasing  the  HONO  group  created  in  6  or  7 

above. 
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9.  The  intramolecule  of  a  hydrogen  from  a  methyl  group  to  the  nitro  group 
breaking  a  C-H  bond  and  creating  an  OH  bond. 

10.  Breaking  the  N-N  bond  and  releasing  the  HONO  group  created  in  9  above. 
14.  Rearranging  the  nitro  group  to  a  nitrite. 

12.  Intramolecular  insertion  of  an  oxygen  atom  from  the  nitro  group  to  a 
methyl  group. 

13.  The  insertion  of  an  oxygen  atom  into  the  methyl  group. 

1A.  The  intermolecular  insertion  of  an  oxygen  atom  and  from  the  nitro 
group  of  one  DMNA  molecule  into  the  methly  group  of  a  second  DMNA  molecule. 
Table  HI  shows  a  summary  of  the  calculations  run  to  date  using  the  MNDO 
and  MINDO/3  methods.  Of  note  is  the  significant  discrepancies  between  the  two 
methods  in  a  number  of  cases.  These  discrepancies  have  not  been  resolved  at 
this  time.  This  work  will  be  continued  under  the  new  work  unit. 
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TABLE  III.  Results  of  DMNA  Calculations 


MNDO 

MINDO/3 

Activation 

Activat ion 

Path 

AH 

Energy 

Ml 

Energy 

1. 

Break  N02 

2.0 

22.5 

23.1 

★ 

2. 

Break  N02(-) 

-21.3 

1.6 

48.7 

★ 

3. 

Break  CH^ 

57.2 

58. 72 

44.8 

* 

4. 

Break  0 

33.7 

94.7 

70.6 

* 

5. 

Break  H 

87.1 

2 

79.1 

87.9 

★ 

6. 

H  Addition 

-63.6 

7.4 

-55.3 

* 

7. 

H  Transfer 

Intermolecular 

21.4 

-v.76 

56.9 

* 

8. 

HONO  Break 

-20.9 

6.5 

25.0 

* 

9. 

H  Transfer 

Intramolecular 

13.0 

<\-66 

24.7 

★ 

10. 

HONO  Break 

58. 91 

* 

14. 81 

* 

11. 

NO2  Rearrange 

-22.4 

68.8 

37.4 

* 

12. 

0  Transfer 

-69.9 

125.4 

-40.8 

* 

Intramolecular 

13. 

0  Addition 

* 

* 

* 

* 

14. 

0  Transfer 

* 

* 

★ 

* 

Intermolecular 

*  Not  Finished 

1.  Product  is  not  stable  and  optimized  geometry  only  approximate 

2.  There  is  probably  no  transition  state  with  these  reactions. 


Chapter  7 

Sulfur  Nitrogen  Calculations 


Introduction 

At  the  present  time  approximately  fifty  molecules  and  ions  (both  positive 

and  negative)  composed  solely  of  sulfur  and  nitrogen  atoms  are  known.  They 

range  in  composition  frm  simple  diatomic  and  triatomic  species  such  as  NS 

2+  - 

and  NNS  to  cyclic  species  such  as  ,  S3N3,  ^2^2  ant* 

extend  on  to  polyhedral  moeties  such  as  S^N^,  S^Ny  S^N^,  etc. 

Ultimately  they  are  all  synthesized  from  the  parent  material,  tetrasulfur 
tetranitride  which  in  turn  is  readily  available  from  the  common  commercial 
chemicals,  sulfur  dichloride  and  ammonia.  Species  are  known  in  all  three 
states  of  matter,  gaseous,  liquid,  and  solid  (including  polymers).  Potential 
applications  exist  in  the  areas  of  liquid  and  solid,  slow  and  rapidly  burning 
propellants,  explosives,  detonators,  electrochemical  battery  applications, 
polymeric  materials  (electrical  conducting  and  as  insulators)  and  in 
materials.  Applications  to  each  of  these  areas  are  outlined  below. 

Thermodynamic  Considerations 

The  literature  has  long  documented  the  explosive  (endothermic)  nature  of 
the  following  materials:  S^N^,  tetrasulfur  tetranitride  (solid)  (thermal 
and  mechanical  instability  in  the  solid  state);  SjNj,  disulfur  dinitride 
(solid)  (thermal  instability);  S^N2>  1-3  tetrasulfur  dinitride,  (liquid) 

S3N3,  trisulfur  trinitride,  and  several  other  materials.  A  recent 
study  at  the  U.  S.  Air  Force,  F.  J.  Seiler  Research  Laboratory  indicated  that 
every  sulfur-nitrogen  compound  known  possesses  a  large  positive  AH^# 
relative  to  molecular  nitrogen  and  sulfur  at  23°C.  Since  AS0  for  most 
decomposition  reactions  can  be  expected  to  be  positive,  reactions  which  lead 
to  a  fragmentation  of  the  parent  molecule,  and  the  production  of  either  1*2, 
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Sg,  or  oxidized  forms  thereof  are  expected  to  be  strongly  exothermic.  The 
observed  explosive  nature  of  these  materials  supports  this  conclusion.  The 
above  listed  materials  offer  a  wide  range  of  potential  solid  state  or  liquid 
state  explosives  or  propellants.  The  mechanical  sensitivity  of 
suggests  application  as  a  potential  detonating  agent.  There  have  been  no 
scientific  studies  to  date  which  attempt  to  evaluate  these  materials  with 
respect  to  the  specific  impulse  available;  interaction  with  specific 
oxidizers,  burn  rate,  etc. 

Electrochemical  Considerations 

A  large  variety  of  positive  and  negative  ions  and  thus  salts  of  sulfur 

nitrogen  compounds  are  known.  They  range  in  structure  from  planar  systems 

2+  4-  +  . 

such  as  ,  S^N,.  and  S^N-j  to  polyhedral  ions 

-  -  2+ 

such  as  ,  S^N  ,  .  They  are  generally  prepared  in 

non  aqueous  solvents  such  as  CH^CN,  DMF,  SbF^  and  others.  Nothing  is 

known  of  their  behavior  in  molten  salt  systems.  There  have  been  a  few 

polargraphic  studies  carried  out.  The  wide  range  of  species  available  suggest 

that  they  might  offer  potential  application  as  supporting  electrolytes.  Many 

of  the  species  possess  ionization  potentials  in  excess  of  10  eV  and  relatively 

low  electron  affinities.  In  the  absence  of  experimental  data  it  is  difficult 

to  be  more  specific,  but  the  diverse  range  of  species  and  structures  available 

suggest  that  this  might  be  a  fruitful  area  of  research. 

Extensive  studies  of  the  electrical  properties  of  polymeric  sulfur 

nitride,  (SN)  have  been  carried  out  and  will  not  be  discussed  here, 
x 

However,  it  should  be  pointed  out  that  synthetic  methods  are  available  for  the 
preparation  of  substituted  materials  of  the  general  type 
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R 


S 


S 


N 


Rl 

N  X 


R*  ,  R  -  CH3,  C2H5,  C6H5 

X  -  1,  2,  3,  4... 


which  can  be  expected  to  show  interesting  properties  such  as  metallic  like 
character,  semiconductivity,  high  melting  points,  etc. 

Application  to  the  Materials  Area 

The  ability  synthetically  to  produce  addition  and  substitution  products  of 
most  of  the  sulfur  nitrides  gives  access  to  a  wide  area  of  materials  of 
possible  diverse  application.  The  area  is  spanned  by  polymeric  materials  such 
as  (HOSN)^  and  its  derivatives,  halogen  substituted  compounds  such  as  NSF, 

^4 ^4 *4*  ^4^3^  an<*  others,  as  well  as  interesting  non 

classical  materials  such  as  Hg5(SN)g,  etc.  Applications  to  the  areas 
oftoxic  materials,  herbicides,  fungicides,  slimicides,  oxidation  inhibitors, 
surfactants,  etc  almost  certainly  exist. 

The  calculations  done  in  this  work  unit  were  in  support  of  those  possible 
applications. 

The  mechanism  of  polyemation  of  N2  was  studied  via  MNPO  type 
molecular  orbital  calculations.  ^  It  was  postulated  that  the  reaction 
proceeds  via  a  cleavage  of  the  S2N2  ring  with  the  formation  of  a  change 
type  intermediate  in  within  a  triplet  or  singlet  state.  Structures  of  and 
activation  energies  were  found  for  the  two  intermediates  and  the  triplet  state 
path  was  found  to  be  energetically  pfeferable. 
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The  isomerization  reaction  HNSO  ♦  HOSN  was  studied  by  MNDO  molecular 
orbital  calculations.  ^  The  S-thiazyl  hydroxide  form  was  found  to  be  lower 
in  energy  by  7  kcal  mol  The  transition  state  was  characterized  with  an 
activation  energy  of  75  kcal  mol  The  electronic  structure  of  the 
reactant,  product  and  transition  state  were  analyzed  using  LMO's. 

The  electronic  structure  of  the  S^N^  anion  was  studied  using 

1  8 

localized  molecular  orbitals  derived  form  G’NDO/2  canonical  M.O.'s. 

Particular  attention  is  given  to  the  nature  of  the  sigma-pi  separation  of 
electronic  charge.  The  calculation  suggested  that  the  anion  is  best  described 
as  a  4*  electron  system  in  contrast  to  the  1 On  model  previously 

2- 

suggested.  The  electronic  structure  of  the  recently  discovered 
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ion  was  studied  by  MNDO  and  LMO  methods.  The  planar  ion  was  demonstrated 

to  conform  to  a  lOir  Huckel-type  moiety.  The  sulfur  trinitride 

bis(triphenylphosphite)  nickel  (II)  dihydrate  compound  was  synthesized  and 

20 

characterized  by  spectral  methods.  It  was  shown  to  be  the  first  example 

of  material  containing  the  SN^  ligand. 

The  general  area  of  the  electronic  structure  of  binary  sulfur-nitrogen 
molecules  and  ions  was  studied  using  density  functional  theory  and  MNDO  types 
of  approximate  wave  functions.  A  detailed  discussion  of  the  identification  of 
equivalent  charge  sites  within  a  molecule  was  presented.  This  work  was 
presented  as  part  of  "Symposium  on  Application  of  Topology  and  Graph  Theory  to 
Chemistry",  University  of  Georgia,  April  1983  and  has  been  published  in  book 
form  as  part  of  the  symposium  volume. 

All  known  binary  compounds  of  sulfur-nitrogen  were  examined  by  MNDO  types  0f 

2  1 

molecular  orbital  calculations  .  The  standard  enthalpy  of  formation  was 

calculated  relative  to  that  of  S.N..  The  materials  were  shown  to  be 

4  4 

highly  endothermic.  Possible  application  to  AFOSR  needs  are  discussed. 
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